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Abstract 

We discuss the existence of asymptotically Euclidean initial data sets to the vacuum 
^ ■ Einstein field equations which would give rise (modulo an existence result for the evolution 

\& ' equations near spatial infinity) to developments with a past and a future null infinity of 

. different smoothness. For simplicity, the analysis is restricted to the class of conformally 

flat, axially symmetric initial data sets. It is shown how the free parameters in the second 
fundamental form of the data can be used to satisfy certain obstructions to the smoothness of 
null infinity. The resulting initial data sets could be interpreted as those of some sort of (non- 
linearly) distorted Schwarzschild black hole. Its developments would be so that they admit 
a peeling future null infinity, but at the same time have a polyhomogeneous (non-peeling) 
past null infinity. 
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1 Introduction 



This article is concerned with providing an example of asymptotically Euclidean, conformally flat 
initial data sets for the Einstein vacuum equations which are expected to have time developments 
for which the two disconnected parts of null infinity have different degrees of smoothness. 

The analysis of the behaviour of the gravitational field in the region of spacetime near spatial 
infinity and null infinity carried out in 15 , 26, 23j — for spacetimes with time reflexion symmetry — 
and in fZ7\ |2E] — for spacetimes without time reflexion — was motivated by the desire of setting 
on a sound footing Penrose 's proposal for the description of the asymptotics of the gravitational 
field — see e.g. [2U1 El EH Penrose's proposal suggests that the asymptotic gravitational field of 
isolated systems should admit a smooth conformal completion at null infinity. 

A large class of spacetimes satisfying Penrose's proposal has been constructed by Chrusciel 
& Delay [S] . The examples provided rely in a refinement of a gluing construction introduced by 
Corvino [5] by means of which it is possible to construct time symmetric initial data sets which 
are essentially arbitrary inside a compact set, but are exactly Schwarzschild in the asymptotic 
region. More recently, a further development in the gluing techniques has allowed to extend 
Corvino's results to the non-time symmetric case, so that initial data sets which are arbitrary 
in a compact region can be glued to a stationary asymptotic region OE]. As a result, all of 
the available examples of spacetimes satisfying Penrose's proposal are stationary in a spacetime 
neighbourhood of spatial infinity. In view of this, it is natural to ask whether initial data sets 
with stationary asymptotic regions are the only ones giving rise to spacetimes with a smooth null 
infinity. In relation to this latter issue, evidence for the following conjecture has been provided 
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Conjecture. The time development of an asymptotically Euclidean initial data set which is 
conformally flat in a neighbourhood of infinity admits a conformal extension to both future and 
past null infinity of class C k , with k a non-negative integer, if and only if the initial data are 
Schwarzschild to order p* where p* = p*(k) is a non-negative integer. If the development admits 
an extension of class C°° ; then the initial data have to be exactly Schwarzschild on B a {i). 

Here, for Schwarzschild up to order p* it is understood that asymptotic expansions of the initial 
data near infinity coincide up to order 0(l/\y\ p ') with those of initial data for the Schwarzschild 
solution — see section 2 for more details. A similar conjecture is expected to hold for non- 
conformally flat initial data — so that the data should be stationary up to order p* if a null 
infinity of class C h is to be attained. 

The calculations leading to the conjecture have shown the possibility of having a spacctimc 
where future null infinity (J^ + ) and past null infinity have different smoothness — note that 

the conjecture requires conditions on both parts of null infinity, in stark contrast with a previous 
version which arose in the analysis of time symmetric situations |26j . The possibility of having 
a spacetime where the two components of null infinity have different degrees of differentiability 
arose for the first time in the post-Newtonian analysis of the relativistic Kepler problem carried 
out by Walker & Will |29l IHH] . Their calculations exhibit a system whose fall-off at null infinity 
is compatible with the Peeling Behaviour at but not at . The latter phenomenon was 
regarded as a hint that requiring smoothness of the gravitational field at past null infinity may be 
a too stringent condition — see e.g. [221 - F° r many physical applications the presence of a non- 
smooth may not be too hampering as long as one can ensure a peeling And even if 
happens to be non-smooth, most of the relevant structure at null infinity can still be recovered 
— see [7]. As pointed out by Penrose: "the issue is not whether regularity at covers all 
situations we would like to call asymptotically flat; it is whether this regularity condition allows 
all the freedom that we need in order to describe isolated systems in general relativity" — see [21] . 

In [22], scripts in the computer algebra system Maple V were used to calculate asymptotic 
expansions of the time development of conformally flat (but non-time symmetric) initial data sets 
in a neighbourhood of spatial infinity. In particular, the components 4>j with j = 0, . . . , 4, of the 
Weyl spinor 4>abcd have an expansion of the form: 

p=|2-j| 

where p is a geodesic distance, and the coefficients <pfi depend on r, an affine parameter of 
conformal geodesies, and on some angular coordinates. In these coordinates r — ±1 corresponds 

to the locus of and J?~, respectively. It turns out that for p = 0, ... ,4, the coefficients 

(v) 

4> 3 extend smoothly through null infinity. However, from p — 5 onwards, one finds a series of 
obstructions to the smoothness of both future and past null infinity. For p — 5 the obstructions 
are time symmetric — that is, the obstructions for J? + and are the same, but at p = 6 
some of the obstructions do not possess this time symmetry. The latter — as it was discussed 
in previous paragraphs — opened the possibility of having a future and a past null infinity with 
different degrees of smoothness. A similar phenomenon occurs, as it is to be expected, if one 
analyses time asymmetric initial data sets which are not conformally flat. However, the overall 
picture is much more involved. 

The calculations in |271 12*5] assume the existence of initial data sets which are expandable 
in powers of 1/r near infinity. Existence proofs for this type of data have been provided in 
|12l 111) . However, the question whether the free parameters in the solution can be chosen such 
that the conditions leading to smoothness-asymmetric null infinities are satisfied remained open. 
This article provides an answer to this question. Ultimately, one would also like to make some 
assertion concerning the existence of the time development. This question is a much more difficult 
one and lies outside the scope of the present work. 

The article is structured as follows. Section 2 is concerned with some aspects of conformally 
flat initial data sets which will be required in our investigations. In particular, some solutions to 
the momentum constraint more general than those of the Bowen-York Ansatz are considered, and 
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Figure 1 : The topology of the initial hypersurface S corresponds to that of time symmetric slices 
of Schwarzschild: two asymptotic regions connected by a "throat" . 

a general existence result for the class of initial data under consideration is recalled. In section 
3 we discuss a certain class of obstructions to the smoothness of null infinity which allow us to 
construct initial data sets whose developments would have prescribed regularity at null infinity. 
Section 4 contains our main result and its proof. 

The present work is a natural extension of the work carried out in |27II28| . We have endeavour 
to follow to the nomenclature and notation of these references, however, for convenience we have 
made use of spin- weighted harmonics instead of certain unitary representations of SU (2, C) when 
performing expansions. Accordingly, our notation has been adapted. The reader in need of a 
remainder on concepts and results on functional analysis (Sobolev spaces, Gateaux derivatives, 
the implicit function theorem) is remitted to [Tl 121 119). 

2 Framework 

We shall restrict our discussion to solutions {h abl Xab) to the vacuum Einstein constraint equations: 



on asymptotically Euclidean 3-manifolds, S, which are maximal — that is, K — — and con- 
formally flat. The metric h ab will be taken to be negative definite. For simplicity, we shall also 
assume that the initial data sets are axially symmetric. Nevertheless, it must be pointed out 
that our analysis can be generalised to the axially symmetric and/or conformally flat settings. 
Furthermore, we shall assume that the 3-manifold S has the topology of the time symmetric slice 
of the Schwarzschild solution: two asymptotically flat regions connected by a "throat — see figure 
1". 

The construction we want to implement depends crucially on the properties of the solutions to 
the constraints in a neighbourhood of infinity. Therefore, we adopt the conformal compactification 
picture. Let S be a smooth, orientable, connected, compact 3-manifold — for definiteness, we can 
think of S 3 . In S we shall consider two special points i and i' — say, the North and South 
poles of S 3 — which will represent the two infinities of the asymptotic ends of the physical data 
(<S, hab, Xab), and so that S = S \ {i, i'}. 

2.1 The conformal method 

Following the standard approach — the Conformal Method — we set 



Let x a denote /i-normal coordinates centred at i. When no confusion arises, we may write x 
instead of x a . In order to ensure the asymptotic Euclidcanity of the initial data we require 



D a xab ~ D b K = 0, 
r-K 2 + K ab K ab = 0, 



(1) 

(2) 



h a b = ^h ab , Xab 



V- 2 i>ab- 



(3) 



V> o6 = 0(\x\- 4 ) as |x| -» 0, 
lim = 1, 



(4a) 
(4b) 
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with \x\ 2 = 8 ao x a x b . Similarly, if y a denote /i-normal coordinates centred at i', we shall require 
that 

^ ab = 0(\y\- 4 ) as \y\ -> 0, (5a) 

lim \y\tf = d, (5b) 
\v\->o 

where d is a positive constant. 

Because of the conformal flatness of the data, we set in a neighbourhood of i, h a b = — Sab, so 
that near i the constraint equations reduce to: 

d a r b = 0, (6a) 
At? = -^ ab r b #- 7 , (6b) 

where A denotes the flat Laplacian. 



2.2 Solutions to the momentum constraint 

The solutions to the momentum constraint in Euclidean space have been extensively studied in 
the literature — see e.g. ^JEJEZ]. For convenience, we will use the version given in [27j, in 
which the solutions to the momentum constraint are derived using the space-spinor formalism 
— see For this, we consider the vector field e%s which in Cartesian coordinates is given by 

x a /|x|, and complete it to an orthonormal basis {eL e%y e ?$\}- Let ip(i)(j) — ipabe^e^^ denote 
the components of the second fundamental form ip a b with respect to this basis. In a space-spinor 
formalism, and under the assumption of a maximal slice, the tensor ip a b will be represented by 
means of a totally symmetric spinor ^Pabcd which is related to V'(i)O') y i a the spatial Infeld-van 
der Waerden symbols vff'- 

= ^f^ijf^ABCD- (7) 

The spinor ipABCD being totally symmetric can be decomposed as 



IpABCD = i>Qe°ABCD + ^ABCD + ^ABCD + ^ABCD + ^ABCD, (8) 

where the spinors d ABCD , i = 0, . . . , 4 are also totally symmetric and 



(l (ABCn)< = <5?( 7.) , (9) 



'4 

--(ABCD), ~ u j U 

the subindex j = 0, . . . , 4 in iabCd) indicating that j of the indices have to be set equal to 1. 
The particular solution to the momentum constraint that we are going to consider can be split 
in six independent pieces: 

IpABCD = i>ABCD + ^ABCD + ^ABCD + ^ABCD + ^ABCD + ^ABCD- ( 10 ) 

The term ^abcd associated with the expansion conformal Killing vector — it gives rise, for 
example, to time asymmetric slices in Schwarzschild 0] EH EH ■ The term ^Pabcd enc °des the 
angular momentum content of the data. While the remaining terms are diverse multipolar terms 
— the first superindex denotes the decay of the term (|cc| -2 , \x\~ x or \x\°) and the second its 
multipolar nature: quadrupolar (2 2 ) or octupolar (2 3 ). For further discussion see |26j . 
In the case of the iPabcd t erm : one has a single non- vanishing component: 



^ = -Oa- ( n ) 



A 



J2 ~ bis 



Assuming axial symmetry, the non-vanishing components of Wabcd are gi ven by: 

t 12 fn 7 12 pK 
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For the remaining terms one has that: 



^ = (-i) J ' +6 - 29 yj 2^J L i,n;« 2-, W, (13) 

for j = 0, . . . , 4. The coefficients £o,n;(j and L^ n . q are freely specifiable complex numbers satisfying 
the reality condition 

Lo,n;q = ^4,n;q- (14) 

The remaining coefficients are calculated via the formulae: 

T {{Li. n , q - L Q , n , 2q )(q+ l)g + 4(n + 7) 2 L ^ q ) , — 

Ll ' n » = (n + 7)(2(n + 8)(n + 6)-( g + 2)( g -l)) ^ + 2 ^~^ ^ 

L 2 , n , q = 3( ^-:g + £4 -"; g) V(g + 2)(g + l)g(g - 1), (15b) 

2(n + 7)^ - g(g + 1) 

T (jMrw ~ Li,n;2 q ){q + l)g + 4(n + 7) 2 L 4 ,» ;g ) /y— ^7 

= (n + 7)(2(n + 8)(n + 6)-(g + 2)( g -l)) V^K^)- (15c) 

In addition, the following regularity conditions will be required to hold: 

^0,-2:2 + ^4,-2:2 = 0, Lo,-l;3 + -^4,-1:3 = 0. (16) 

The latter, together with (|14J) that Lo,-2;2> -^4,-2;2) io.-i;3 and £4,-153 are pure imaginary. 
The conditions (|16|) are intended to preclude the appearance of a certain type of logarithmic 
divergences at the sets where null infinity touches spatial infinity — see |151 177| , and also the 
discussion in section 3. 

From the point of view of the normal coordinates y a centred on i', the different terms of 
4>abcd have an analogous decay to the one near i: 

^BCD = 0(\y\~ 3 ), r ABCD =0{\y\-% rA k BC D=0{\y\ n ). (17) 

Thus, all the terms — save for 4> < abcd — are singular at both i and i' — see |12| . 
For latter use, it is noted that if ipab is of the form given by equation I jlOjl then: 

Phi 




<■„,,<•"'■= •:\ucdc au '' u = — [ '— ■ 12 J- -!- -l\ K I^.J 2 Y m \ • (18) 

Note that the above expansion truncates at order C(l/|a;| 2 ), the first contributions of the higher 
multipolar terms in ^™£ fe appearing at order C(l/|a;| 5 ). Also, tp a bilJ ab < 0. Finally, we note that 
for a neighbourhood B a {i) of radius a centred at i, one has that 

\x\ s ip ab r b eE™(S\ {»'}), (19) 

where E°°(B) = {/+ \x\g \ f, g £ C°°{B)}. Similar behaviour is obtained on for \y\ 8r ip a bip ab ■ 

2.3 General existence results 

The assumption on the conformal flatness of the 3-metric h a b and the choice of solutions to the 
momentum constraint made in the previous section fulfil the hypothesis of the existence results 
obtained in Thus, one has: 

Theorem 1 (Dain &: Friedrich, 2001). Ifh a b is conformally flat, and\ab is given via equation 
U(fy . then there is a unique, positive solution $ to the Licnerowicz equation Ifibp with the boundary 
conditions and \5b\l . Moreover, by virtue of )iy\) one has that on S \ {«'}.' 

= n + w, (20) 

\x\ 

where W £ E°°(S\{i'}). 
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Because of W S E°°(S \ {«'}), the function W can be expanded in a neighbourhood of i as: 

W = j + (d a x a + m'\x\) + (q ab x a x b + d' a x a \x\) + (o abc x a x b x c + q' ab x a x b \x\) + 0{\x\ A ), (21) 

where m is the ADM mass of the initial data, and m', d a and d' a , q a b and q' ab , o abc are, respectively, 
monopolar, dipolar, quadrupolar and octupolar terms. In reference |2(j) it has been shown that, 
without loss of generality, one can make use of the translational freedom in the setting to remove 
the dipolar term, d a , in the previous expansion. Proceeding in such a way, we are left with 

W = j + m'\x\ + (q ab x a x b + d' a x a \x\) + (o abc x a x b x c + q' ab x a x b \x\) + 0(\x\ 4 ). (22) 

The explicit calculations carried out in show that in fact m! = and d' a = 0. 

Alternatively, one could have used the spherical coordinates associated with the normal coor- 
dinates x a . Due to our assumption on axial symmetry the expansions are given by 1 : 

w = y + h r " £ ^ lYw + aV 3 £ ^ lYl ° + ° (r4) ' (23) 

1=0 ' 1=0 

where r = \x\, and LUp-i_ m S M. 

The function W satisfies the elliptic equation: 

AW = -\^ ab r\ll\x\ + W)- 7 . (24) 



Thus, the coefficients 0J2-i. m and 0J3-^ m are determined by the boundary conditions I|4a|l - (j5b|l 
and the parameters of ip ab — that is, A, J , io,n;fe and L^^-k- Thus, sometimes we shall write 
W = W[A, J, L](x). The elliptic nature of (|24|l reveals that W contains information of global 
nature. The local part — around infinity — is encoded in the l/\x\ term of (|20[1 corresponding to 
the conformal flatness of the data. 

From equation (|2~4fl and taking into account the boundary conditions lf4"af> - lj4b() one gets the 
integral representation: 

j_ r Mx>)r b {x>) i 

4tt y R 3 8(l/|x'|+W(x')) 7 \x-x'\ ' 1 ' 

The function — l/47r|x — x'\ corresponds to the Green's function of the Laplacian on IR 3 . 

For later use, we show that the function W = W[A, J, L](x) depends — at least — in a C 1 
fashion on the parameters A, J £o.n;fc and L^^-k which determine the tensor ip ab . 

Lemma 1. The solution W of equation \2J$ , with ip ab given by equation $1U\) is C' 1 in the 

parameters A, J £ R and Lo,ri;fc, E^.n-.k G C 

To see this, we rewrite equation (|24|) as 

G(J; W) = AW + lip ab ^ ab (l/\x\ + W)- 7 = 0, (26) 

o 

and concentrate on the dependence of W with respect to J — the analysis with respect to the 
other parameters of ip ab is similar. The mapping Q : Rx W 2,P (S) — ► L P (S), p > 2, where W 2 ' P (S) 
denotes the Sobolev space of functions on 5 with L p -integrable weak second order derivatives, is 

reference 1271 . on which the present work is based, a certain type of unitary representations of 5(7(2, C), 
T m j, is used. Here, instead, we make use of the more standard spin- weighted spherical harmonics s Y nm . The 
relation between the two sets of functions is given by 



V . f!\s+2n—m / 2n + 1 ^ n — m 



n — s 



Accordingly, we have slightly changed our notation with respect to that of 1271 . 



(i 



at least a C 2 mapping (in the LP norm) between Banach spaces — actually, it can be shown to be 
C°° , but C 2 will be enough for our applications. That Q € C°(M x W 2,p , L p ) can be shown using 

\\G(J,W)-g{Jo,W )\\ L P < \\AW- AWo\\lp 

+C\\(1 + \x\W)- 7 \\ LP \\\xU ab r b (J) - \x\ 7 ^ ab r b (Jo)\\ LP 

+q|| a; r^ ah (Jo)l|L P ||(i + \x\wy 7 - (i + \x\w )- 7 \\ LP , (27) 

and noting that AW converges to AWq is W goes to Wo as W 2 ' p is a Banach space, and that 
ipab'4' ab 1S m our case analytic function of J. To establish the C 1 character of the map Q one has to 
perform a similar discussion with the Gateaux derivative of G, dcG '■ x W 2 ' p x W 2,p — > L 2 ' p , 
which is given by 

d G G(J,j; W, w) = Aw- \jdj^ ab r h {l/\x\ + W)~ 7 + 7 -w^ ab {l/\x\ + W)~\ (28) 

Again, the analyticity of ip a b with respect to J is crucial to establish the continuity of &gG- 
Higher order derivatives are dealt with in a similar fashion. 

Now, let us consider a solution Wq to equation (|21f> . corresponding to a certain value, Jo, of 
the angular momentum, so that 

G(Jo,W ) = 0. (29) 
The existence of such a solution is guaranteed by the theorem 1. The linearised operator 

Cw = Aw - l^ ab (l/\x\ + W)~ 8 w, (30) 

o 

£ : W 2,p — > L p is an isomorphism: due to ipab' l P ah > the only solution to Cw = is w = (seen 
using integration by parts); further, C is selfadjoint and, thus, using the Fredholm alternative 
it is also surjective — see eg. |12| . So, from the implicit function theorem we know there exists 
solution W — W[J]{x) to equation 126|) for J sufficiently close to Jo — cfr. [I]. By unicity, this 
solution has to coincide with that one given in theorem 1. More importantly for our purposes 
is that W € C 2 {U,X) —actually C*°°!— for some U C R, Jo £ U and X C W 2 > p , W(J ) G X. 
Poinwise differentiability with respect to J follows from the Sobolev embedding theorem: if 
(k - r - a)/n > 1/p then W k ' p (U) C C r ' a (U), with U C M™ —see, for example [2- In our 
case k = 2 (at least!), n = 1, r > 1 and p suitably large — say, p > 2. Note that if one only 
had W £ C 1 (t/, X) then the pointwise continuity — with respect to J — of djW can not be 
concluded. 

3 Obstructions to the smoothness of null infinity 

As mentioned in the introduction, the calculations carried out in |27j provide asymptotic expan- 
sions for the components of the Weyl spinor which allow to relate the structure of the initial data 
with the radiative properties of the development. Here, and in what follows we will concentrate 
our attention on the null infinity associated with the asymptotic end i. What happens at the 
null infinity assigned to the second asymptotic end, i', will not be of relevance for our purposes. 
Assuming axial symmetry, the components of the Weyl tensor can be written as: 

& = E E a ^(r)r p 2-jY l0 , (31) 

P=0l=\2-j\ 

where j = 0, . . . ,4. The coefficients ckj jP; ;(r) for p = 0, . . . ,4 are polynomial in r. 

For p = 5, they are polynomial in r for I = 0,1, 3. However, for I = 2 one has that 

a„ 5;2 (T) = T 5;2 ((1 - T) 7 - j Vj{r) ln(l - r) + (1 + r) 3 +^ 4 -,(r) ln(l + r)) + Q(r), (32) 

where Vj(r), V^-j(r) and Q(t) are some polynomials in r. In particular, Vj(±l) ^ 0, p4_j(±l) 7^ 
0. 
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Similarly, for p = 6 the coefficients ctj lP - t i(r) are polynomial for p = 0, 1, 4, while the coefficients 
of the / = 2 harmonics are of the form 



«i,6 5 2 = T+ 2 (l - r) 8 "' ln(l - t)V 3 (t) + T^ 2 (l + r)^P H (r) ln(l + r) + Q(r), (33) 

where again Vj(r), Vi-j{r) and Q(r) are some polynomials in r, P,(±l) 7^ 0, ^-^(±1) 7^ 
— different to those in equation l|32[l . Finally, 

a jA3 (r) = T 6;3 ((f-r) 8 -^ 3+ i(r)ln(l-r) + (f + t) 4+j '7V,(t) ln(l + r)) + Q(r). (34) 

The coefficients aj jP; g(r) for p > 7 are expected to exhibit a similar pattern. In particular, it 
is conjectured that there is an infinite hierarchy of coefficients associated with logarithmic 
terms of the form discussed above. These details will not be relevant for our analysis. 

The coefficients Ts ;2 , T^ 2 and Tg ; 3 — the obstructions to the smoothness of null infinity — 
are given in terms of the freely specifiable data by: 



T 5;2 = 9\ 

L 6;2 



1 6;2 — 



3 2 37602 „ 3099 r- , 2046 r- . n 448 ,-. r 

9y -m ^ 2;2i2 + -^-mj - ^"V6m ia + — Vemift + ^ V6ifr, (35) 

2198208 o , 62451 , ~ 62691 /-. .„ 1791 r- n , 263327.,, 

^M.-'-f^^V-m,, (36) 

„ ff 565753248. 36399 r n 3808 r r 7272 ^- , . 

- lV^>.» - -82585-^ + HT^' + 415"^^ " 415 ^ (3?) 

where in order to render the above expressions more readable we have set 

^0,-2:4 = 1", ^4,-2:4 = — ia, (38a) 

A),-i;4 = Pr + Li,-i;i=0R-i0i, (38b) 

^0,-1,6 = i7> i4,-l;6 = -i7) ( 38c ) 
^0,0;4 = #R + i<5/, ^4,0,4 = &R - i^J, (38d) 

in accordance with the reality conditions l|14|l and the regularity conditions l|16|) . 

The role of the obstructions l|35|) - (|37|l can be better understood by noting that the conformal 
factor — rendering the representation of the region of spacetime near null and spatial infinity 
in which the components of the Weyl tensor l|31|l are calculated — is of the form, 

where D a dD a d /i? 3 = 0{p). The locus of null infinity is given by J ± = {t = ±l,r ^ 0}. In 
addition, we define I = {p = 0, |r| < 1} and 1^ = {p = 0, r = ±1} corresponding, respectively, to 
the cylinder at spatial infinity and the critical sets where null infinity "touches" spatial infinity. 
Thus, the coefficients ctj^qir) in the expansion (|31|l . considered as functions on / diverge on the 
critical sets 1^ if the relevant obstructions are not satisfied. Because of the hyperbolic nature of 
the equations governing the evolution of the gravitational field, one expects that these singularities 
will propagate along null infinity. 



3.1 A particular choice of the data 

Inspired by the above discussion, we are interested in constructing an initial data such that: 

T 5;2 - T+ 2 = T 6;3 - 0, (40a) 
T 6;2 ^ 0. (40b) 

The evolution of such an initial data set will possess — modulo an existence proof for the conformal 
field equations which is valid up to the critical sets, I , which is not yet available, see e.g. |17j— 
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a null infinity where J'^ and J?~ have different smoothness. Assuming that the conditions 140a(l 
and (|40b|> are sharp, a rough count using ll3*2"|l - (f3~4"|) suggests that the time development of our 
initial data will have an which, at best, will be of class C 3 and an of at least class C 4 . 
In particular it is claimed the following: 

Conjecture. The time development of an initial data set (h a b,Xab) such that the conditions \4-0aty 
and hold admits a peeling future null infinity and a polyhomogeneous past null infinity. More 

precisely, using and an adapted gauge the components of the Weyl tensor 2 decay, near <y + , as 

fa = 0(l/r), fa =0{l/r 2 ), fa = 0(l/r 3 ), fa = 0(l/r 4 ), fa = 0(l/r 5 ). (41) 

While near the decay will be 

fa = 0(l/r), fa=(D(l/r 2 ), fa=(D(l/r 3 ), fa = 0(l/r 4 ), fa = 0(lnr/r 5 ). (42) 

That is, the resulting spacetime is expected to peel at J? + , but not at . For a discussion 
of the notion of polyhomogeneity at null infinity and non-peeling spacetimes see for example . 

4 The main result 

In this section it is shown how the freely specifiable parameters contained in the initial data 
discussed in section 2 — the scalars A, J and a, j3 R , j3j, 7, 5 R and Sj — can be used to construct 
initial data satisfying the conditions (|40a|l and (|40bjl . To this end, one has to find a way of 
controlling the coefficients lu 2 -2 and 103-3 hi the expansion (|23|l . 

Theorem 2. Given Aq € K, Ao 7^ and a conformally flat initial data set with second fun- 
damental form given by 1) lUp and satisfying the boundary conditions \4b^ and f5b\) . there ex- 
ists a neighbourhood 13 C R 3 , (0,0,0) € 13 such that if i^J , (3 R , (5i) € 13, then there are unique 
A = AiJ^R^j), a = a{,J,p R ,Pi), S R = 6 R (J, p R , (3i), 5j = 5j(J, [3 R , fc), 7 = i{J,p Ri fa) 
with A{Q, 0, 0) = Ao, a(0, 0, 0) = S R (0, 0, 0) = 5/(0, 0, 0) = 7(0, 0, 0) = 0, such that 

1 5:2 — 1 6;2 — 1 6;3 — U, 

The condition 

Te ;2 * 0, 

is satisfied ifJ^Q,l3 R ^ 0, /?/ ^ 0. 

Thus, given a particular time asymmetric, conformally flat Schwarzschild initial data set, 
there is always a (non-linear) perturbation of it — actually an open set of them — satisfying the 
conditions (|40a|l and l|40b(l . Note that these perturbations have necessarily non- vanishing angular 
momentum; furthermore, from a careful analysis of the proof it stems that all the parameters of 
the second fundamental form have to be non-zero — justifying the somewhat complicated-looking 
choice of the extrinsic curvature made in 2.2. The different parts of the proof will be discussed 
in the sequel. 

4.1 Integral representations of lo 2 -2 and a> 3;3 . 

Our first task is to obtain some expressions which enable us to control the expansion coefficients 
in The explicit calculations of ,27, show that: 

1 1 3 

W = - + ^r 2 cj 2l2 Y 20 + -r 3 £ cj 3;l Y lm + 0(r 4 ), (43) 

' 1=0 

2 The convention used here to denote the different components of the Weyl tensor is exactly the opposite of the 

standard Newman-Penrose (NP) convention. That is, our tf>o corresponds to NP's 1/14, and so on. 
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in S \ {i'}, where 

/3J 2 A 2 \ 

UW = -2Vt[- T + ggj, (44a) 

W3ji = J 2 , (44b) 
oj 3 . 2 = 0. (44c) 

On the other hand, the coefhcients w 2;2 and W3 ; 3 do not admit such closed and neat expressions. 
They contain information of global nature. Thus, in order to control them we have to resort to 
the integral representation (|25|) . 
We begin by noting that: 

d rr W = w 2;2 y 20 + 0(r). (45) 

In order to differentiate the integral representation l|25|> we have to proceed with care, for the 
Green function l/\x — x'\ is singular if x = x' . We begin by letting 



F(x) = (46) 



Due to the existence results of Theorem 1, F(x) G E°° (S \ {«'}), and F(x) — 0(\x\). Accordingly, 
d^F, k = 1, 2, 3, are regular at r = \x\ = 0. And so, 



w = w- L^r^-r^*'' 

327tJk 3 \x-x'\ 



32tt ,/lR 3 

Frome where 



di k) W = [ s d^F(x - x')-±-d 3 x', 
32tt Jk 3 \x'\ 



L F (x-x')Ad S x'. (47) 



(48) 



^2;2 = / (Y 20 \d rr F{-x'))—d 3 x', (49) 

327T./K 3 x' 



(Y lm \H(0,ip)) = / Y lm (e,<p)H(e,<p)em0d0d(p, (50) 
Jo Jo 



with fc = 1, 2, 3. Consequently, we find that 

i r „. ,„ „, , i 



where 



denotes the component of H on the (Z, m)-harmonic. Similarly, one finds that 

^3;3 = ^T I (Y 30 \d rrr F(~x')}-±-d 3 x'. 

327TJK 3 \x'\ 



(51) 



4.2 Satisfying the conditions for a "smoothness asymmetric null infin- 
ity". 

In order to see how the conditions 14()a(l and l|40b|l can be satisfied, we begin by considering 
imaginary parts of the obstructions 1|35JI - 137JI : 

Im(T 5;2 ) = -^V6m 2 a + ^^mfr + gjj^ = 0, (52a) 

. 62691 263327^ 116 . r 183184 „ . . 

Im T+ = VQABj H 3i V6.46/ V6<5j = 0, 52b 

V 6 ' 2 > 2408 M 31605 M 43 13545 V ; 

_ . 565753248 , 7272 , . 

Im T 6 . 3 = .J 3 V307 = 0. (52c 

V b ' iJ 82585 415 ' y ' 
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The condition (|52a[) can be used to solve for a in terms of (3i and 5[. Condition l|52b|l can be 
used to solve for Si in terms of A and /?/. Finally, condition (|52c|l can be used to write 7 in terms 
of J, rendering 7 oc J 3 . On similar lines 

„ ,„+ \ 2198208 - 62451 1791 f7.n 1791 

Re(T+ 2 ) = 6Q65 J 2 + ~^A7 2 + -j^VQPr - —VQ5r = 0, (53) 

can be used to solve for 8r in terms of A, J and (3 R . Hence, we are left with only 4 free real 
parameters: A, J, (3r, 

Satisfying Re(T5 ;2 ) = and Re(T6 ; 3) = is much more subtle. For this we turn to the 
integral representations l|49|l and 151|) of the coefficients w 2; 2 and tu^ of W. Consequently, one 
finds that 

Re(T 5;2 ) = ^-^m 2 j (Y 2Q \d rr F(-x'))±-d 3 x' + ^mj 2 = 0, (54a) 
AZlT V 7T 7jj 3 \x'\ 199 

Re(T 6;3 ) = f (Y 30 \d rrr F(-x')) -±-d 3 x' + ^VGJPi + ^VGJSj = 0. (54b) 

on V 7r JM. F I 415 415 

Now, we note that W = W[A, J, (3r,Pi](x). That is, the function W depends on the the choice 
of the free parameters A, J , (3r, (5j. For example, given our choice of boundary conditions, if 
A = J = (3r = [3i = 0, then the initial set reduces to that of the standard time symmetric, 
conformally flat slice of the Schwarzschild spacetime, and accordingly W = m/2. More crucially 
for our purposes, if A = Aq 7^ 0, but we keep J = (5r = Pi = 0, then the initial data corresponds 
to a conformally flat, time asymmetric slice in the Schwarzschild spacetime, and — see |27|— 

W=^-^A 2 \x\ 2 + <D(\ X n (55) 

The latter fact suggests the use of an "implicit function theorem" -type of argument to see if there 
are some values of A, J , (3r and /3/ such that (|54a|l and (|54b|) hold. 

Given that the choice A = Ao 7^ 0, J = (3r ; = (3j = corresponds to Schwarzschild data, we 
have that 

Re(T 5;2 )[.4 , 0,0,0] = 0, (56a) 
Re(T 6;3 )[A, 0,0,0] = 0; (56b) 

the smoothness of of static spacetimes at the critical sets where null infinity "touches" spatial 
infinity has been discussed in JSj. On the other hand, using lemma 1, Re(T 5;2 ) and Re(T 6; 3) 
depend in, at least, C 1 fashion on A. From l|18[) and (|55fl it follows that 

^Re(T 5;2 )[A,0,0,0] ^0, (57a) 

cl4Re(r 6;3 )[A>, 0,0,0] 7^0. (57b) 

Thus, if we set 

F(A,J,p R ,/3 I ) = (Re(r 5]2 )[A,J,/3 R ,p I }^ + (Re(T 6 , 3 )[A, J,p R ,0if) , (58) 
one has that 

^"(A), 0,0,0) = 0, (59a) 

d A F(A , 0,0, 0)^0, (59b) 

so that from the implicit function theorem — see e.g. ^ HOj — there exists a neighbourhood 
13 C R 3 , (0, 0, 0) € 15 and a C 1 function A{J, /3 R , /?/), .4(0, 0, 0) = A such that 

F(A(J,/3 R ,/3 I ),J,p R ,f3 I )=0, (J,0 R ,0i) £15. (60) 

The latter proves essentially proves our main result. A careful inspection reveals that T^ 2 7^ if 
J 7^ 0, /3r 7^ and /3j 7^ 0. We note that a similar argument with the choice A = J — (3r — (3i — 
does not work for it implies dx^ab^ ah — d j^W = so that c^-^O, 0, 0, 0) = — the implicit 
function theorem can not be applied. That is, the Schwarzschild initial data that one is perturbing 
have to be time necessarily asymmetric. 
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5 Conclusions 



The initial data we have constructed in this article is a non-linear perturbation of time asymmetric 
Schwarzschild initial data. It is conformally flat and has non-vanishing angular momentum. 
Its time development should give rise to a radiative spacetime which could be interpreted as a 
distorted Schwarzschild black hole. The crucial remaining problem is to obtain an existence proof 
for the evolution equations which is valid up to (and beyond) the critical sets where null infinity 
"touches" spatial infinity — see the discussions in (El Q2| • Obtaining such a proof will require 
deeper insights into the structure of spatial infinity than the ones currently available. 
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